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Abstract A fast-wave pulse in a simple, cold, inhomogeneous MHD model plasma 
is constructed by Fourier superposition over frequency of harmonic waves that are 
singular at their respective Alfven resonances. The pulse partially reflects before 
reaching the resonance layer, but also partially tunnels through to it to mode convert 
to an Alfven wave. The exact absorption/conversion coefficient for the pulse is shown 
to be given precisely by a function of transverse wavenumber tabulated in Paper I 
of this sequence, and to be independent of frequency and pulse width. 
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1. Introduction 



Recently, Cally and Andries (2010 I (Paper I) explored mode conversion and resonant 
absorption in a simple model consisting of a uniform magnetic field oriented in the 
2-direction threading a cold magnetohydrodynamic (MHD) plasma with a monoton- 
ically decreasing density stratification in the .i-direction. The Alfven speed a(x) is 
consequently monotonically increasing in x. A monochromatic, plane, fast wave with 
frequency uj and y and z wavenumbers k y and k z incident from x — > — oo classically 
reflects where uj 2 = a 2 {k 2 + k 2 ), but in fact partially tunnels through to the Alfven 
resonance uj 2 = a 2 k 2 located distance X to the right of the reflection point. We adopt 
the simple exponential Alfven profile a = a u e x ' 2ft , which may be regarded as deriving 
from a local tangent approximation to In a in the neighbourhood of the resonance. 
It was shown for this profile in Paper I that the fraction A of the incident wave- 
energy flux that is absorbed by the resonance and converted into an Alfven wave 
depends solely on the two parameters 9 = tan~ 1 (k y /k z ) and a = (k ± h) 2 / 3 sin 2 9, 
where k± = \Jk 2 + k 2 . The remainder of the energy is reflected back to the left, 
also in the form of a fast MHD wave. Without loss of generality, the density scale 
distance h is set to unity throughout. 

Intriguingly, absorption may not exceed A Q = 49.37% (see Figure [T}, the max- 
imum corresponding to a = a = 0.464434 and 9 — > 0°, and falling rapidly as a 
increases or decreases from this value. Of course, k ± — )• oo in the limit where 9 
vanishes but a remains finite. As is well known, the Alfven resonance corresponds to 
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Figure 1. Absorption coefficient A(cr,8) for a? <x e x / h and 9 = 5°, 10°, . . . , 85° (coloured 
curves, top to bottom). The heavy dashed curve corresponds to 8 = 0°. The heavy vertical 
line is at a = ctq = 0.464434, where SA(o, 0) attains the maximum value Aq = 0.493698. J? is 
tabulated in Appendix C of Paper I. [This is a reprinting of Figure 2a of Paper I.] 



a singular point of the governing differential equation, where the wave-energy density 
is not only unbounded, but also not integrable. This is because energy is assumed 
to have been absorbed at constant rate over an infinite time. 

The existence of the resonance is a consequence of the magnetic field being perpen- 
dicular to the inhomogeneity direction x (assuming specified Fourier dependencies 
in the y and z directions). We choose this special case precisely because we wish to 
explore the behaviour of a wave pulse meeting a resonance in the simplest possible 
model. With oblique field, mode conversion still operates and may exceed Aq. In the 
solar context, the case we are considering is most relevant to horizontal magnetic 
field in a vertically stratified atmosphere. 

To better explore the nature of the resonance, and the relationship between 
resonant absorption and mode conversion, Cally and Andries (20101 constructed 
a well-behaved spatial ray bundle from the singular wave solutions by Fourier com- 
position. The modes from which this bundle was synthesized all had the same u 
and k y but were distributed in k z as a gaussian about a central k zQ . This had the 
dual effects of both localizing the conversion/absorption region in z and smearing 
the resonance. A very clear mode conversion process became apparent graphically, 
with the conversion coefficient being an appropriate average of the ^ values drawn 
from Figure [T] 

On presenting these results in seminars at several institutes, one of the authors 
(PSC) was consistently asked about how this would generalize to temporal wave 
packets; i.e., a sharp (or otherwise) incident wave burst with single k y and k z . This 
is relevant to waves excited by discrete events in the solar corona for example. The 
answer is quite simple, and in the interests of completeness is presented here. For 
conciseness, necessary results from Cally and Andries ( 2010 \ will be assumed. Further 
discussion of context and of the related literature may also be found in Paper I. 

Whereas the spatial ray bundle case deposits energy continuously into the res- 
onant layer in the vicinity of a localized height [20] and mode conversion can be 
identified as an Alfvenic-wave energy flux on z > zq taking energy away from this 
site, the temporal ray packet instead deposits energy equally at all z along the 
resonance layer but only for a short time. The issue is then the total energy per unit 
y— z area left on the resonance once the driving fast wave has departed, as a fraction 
of the integrated incident flux. We ask how that fraction depends on frequency, 
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wavenumbers, and packet duration, and in particular if it can in any way exceed 
A = 49.37%. 



2. Temporal Wave Packet 

Let x{x,u) be the (singular) monochromatic solution for the dilatation \ = V>£ 
corresponding to frequency u, where £ is the plasma displacement. Assuming a x = 
Xexp[i(k y y + k z z — uit)] transverse and time dependence, the governing equation is 
(c./., Equation (3) of Paper I) 



d 

da; 



1 d\ I / ., k 2 



oj 2 /a 2 k 2 — 1 dx 



+ t^ 2 -^%nJx = o. a) 



Note that x has a logarithmic singularity at the Alfven resonance u 2 = a 2 k 



z' 



A temporal ray packet may be constructed from x by Fourier composition, 

1 (°° 

X{x,t) = —J e-^ t A{uj) X {x,uj)Au (2) 



where A(u) = v / 2~7rTexp — \ (oj — u>o) T 2 is a gaussian distribution about uj in 
the frequency domain. T is the corresponding temporal width of the packet. The 
harmonic spatial dependence in x is implicit: x( x , V, z , t) = x{ x > t) exp[i(k y y + k z z)\. 

The analysis of Paper I, in particular involving Equation (15), shows that x 
depends only on a and 6, and not on frequency to except for a trivial translation: 

X(x, w) = Xae {x - In ^ ^ , (3) 

where the reflection point for the central frequency loo has been arbitrarily placed 
at x = 0. The resonance is then at X(u) = lnsec 2 # — ln(w 2 /a>Q), which of course 
is smeared by the superposition of the continuous range of frequency in the packet. 
Unlike the spatial ray bundle case (where X was smeared by the range of 9), a single 
numerical solution x (for given a and 9) suffices to construct x(x,t). 

Since the elemental absorption coefficient J? depends only on a and 9 (and thus 
on k y and k z ), and not on frequency oj, it is clear that the temporal wave packet con- 
structed from waves with a common a and 9 but a range of u must have absorption 
coefficient Si{a, 9) precisely, independent of ui and T. We may therefore definitively 
answer the question of whether absorption may exceed A in the negative without 
further ado. Nevertheless, we illustrate this conclusion with a numerical example. 

Having constructed x> the plasma displacement £ = (f x , £ y , 0) may be formed by 
integrating the driven Alfven equation 

^^f +/c ^ = V pX, ( 4 ) 

where V p is the x-y part of the gradient operator. This clearly shows how V P x acts 
as a source term for Alfven waves, since A = a~ 2 d 2 /dt 2 + k 2 is the free Alfven 
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Figure 2. Left: Energy density E as a function of x and t for the case discussed in the text. 
The fast ray bundle enters the box at x^ = —3 at around t = — 9.2 and leaves at about t = 9.2, 
leaving 46.4% of its energy in the resonant layer. The vertical green line is the resonance for 
the central frequency u)q = 9. Right: Rex (background shading) and Re^y (superposed red 
and blue curves). 



operator. Equation Q has the exact solution (assuming £ vanishes as t — > — oo) 

/t 
V p x(x,t') smak z (t-t')dt' . (5) 
-oo 

We require expressions for the wave-energy flux [F x ] in the x-direction and for 
the wave-energy density [E]: 

F x = Re{- X *v x } , (6) 

where v x = d£ x /dt is the x- velocity, and 

E = \ [(K| 2 + \v y \ 2 )/a 2 + k 2 M x \ 2 + \Q 2 ) + \ X \ 2 ] . (7) 

A factor B 2 /n has been suppressed in both, where B is the background magnetic-field 
strength and /i is the magnetic permeability. The three parts of E are respectively 
the kinetic, the magnetic tension, and the magnetic pressure energies. Assuming that 
the incident wave enters the system x > xl (where xl < 0) entirely during t < 0, 
and leaves during t > 0, the total energy per unit y-z area entering through a left- 
hand boundary at xl is E + = J_ F x (xl, t) dt > 0, whilst that leaving through 
the same boundary after reflection is £L = F x (xl, t) dt < 0. The absorption 
coefficient is therefore (E+ + E-)/E + , or alternatively lim^oo J°° E(x, t) dx/E + . 
The equivalence of these two expressions is a non-trivial check on the numerics and 
is attained to good accuracy in our calculations. 

The case we examine is k y = 3, k z = 9 (and hence a = 0.448, 9 = 18.43°, 
J% = 0.4638, quite close to the maximum possible absorption), T = 1, with uj = 9. 
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A total of 2 9 x 3 = 1536 Fourier modes are employed, and the transform (p| 



performed numerically using the fast Fractional Fourier Transform algorithm ( Bailey 



and Swarztrauber, 1994), implemented in parallel. The energy density and displace- 



ment vectors are plotted in Figure [2] The red curve is the exact ray path predicted 
by classical ray theory for the central frequency ui ; viz. ljq i 2 = Ak^(e~ x — 1). The left 
frame shows clearly that the incident fast wave splits into two near the resonance, 
partially reflecting as a fast wave but also partially remaining on the resonance layer 
as an Alfven wave. We confirm that (E + + EJ)/E + = 0.4638 as expected. The Alfven 
wave is characterized by its transverse displacement [£ y ] explicitly shown in the right 
frame, with \ representing the incident and reflected fast wave, which is compressive. 
The increasing inclination of these £ y wave fronts with time denotes phase mixing. 
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Figure 3. Left: Integrated total energy (full curve), Alfvenic energy (long dashed) and 
non-Alfvenic energy (short dashed) per unit y— z area in the region —3 < x < oo as a function 
of time t. Centre: Alfvenic (long dashed) and non-Alfvenic (short dashed) energy densities as 
functions of x at t = 0. Right: Alfvenic and non-Alfvenic energy densities at t = 15. 

Figure [3ji displays the total energy per unit area in i > —3 showing clearly 
how it rises monotonically as the wave enters the box from the left, stays constant 
for a period, then decreases monotonically as the reflected fast wave departs but 
nevertheless leaves a substantial portion of its energy (the Alfvenic part) on the 
resonance in perpetuity. This, together with Figures [3J5 and [3ji, confirms that the 
Alfvenic energy [Eau = |(|w y | 2 /a 2 + fc 2 |£j,| 2 )] is generated around t = and is 
the only significant energy remaining on the resonance layer once the fast wave has 
departed. 

Finally, Figure [4] and the accompanying animation show the fast wave as a shaded 
avocado rendering of x over-plotted with a red-blue represention of the Alfven wave 
in the form of £ y . They illustrate the reflection and mode conversion very clearly, as 
well as the phase mixing of the Alfven wave. The movie in particular gives a good 
physical impression of how the refracting fast wave lines up with the Alfven wave 
and gives it a "push along" in the conversion region. This coincidence of the phase 
velocities of the converting and converted waves is fundamental to the process. 



3. Conclusions 

We have once again illustrated that constructing a ray packet - temporal on this 
occasion - from singular resonant monochromatic wave solutions provides a clear 
picture of fast-to- Alfven mode conversion in a cold MHD plasma, with the conversion 
region localized in both x and t. The solutions are smooth everywhere. We find that, 
since the absorption/conversion coefficient depends only on transverse wave numbers 
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[ky] and [k z ], and not on frequency [a;], such a bundle constructed out of waves of 
different frequency but the same k y and k z suffers exactly the absorption predicted 
by Figure [T] and tabulated in Paper I. In particular, it is not possible to engineer a 
higher conversion coefficient than Aq = 0.4937. Of course though the frequency is 
crucial in determining the position of the Alfven resonance. 

Although the model adopted here is extremely simple, it is instructive when con- 
sidering impulsive waves in the solar corona. We may surmise that fast MHD waves 
generated by flares and other such localized events may easily end up partially mode 
converting to Alfven waves in resonant layers as they propagate into stronger-field 
or lower-density regions, and that the width of the pulse should have little influence 
on the ultimate absorption/conversion fraction, which may not exceed about 50% 
when the magnetic field is oriented perpendicular to the direction of inhomogeneity. 
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